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Physical  and  Optical  Properties  of  Falling  Snow 

GARY  KOH 


INTRODUCTION 

Many  military  electro-optical  devices  detect 
visible  and  infrared  radiation  that  is  emitted  or  re¬ 
flected  by  distant  targets.  The  performance  of  these 
electro-optical  devices  that  rely  on  the  propaga¬ 
tion  of  optical  energy  between  two  points  in  space 
is  largely  dependent  on  the  atmospheric  condi¬ 
tions.  Snow  in  the  atmosphere  scatters  and  ab¬ 
sorbs  visible  and  infrared  radiation.  1  his  reduces 
the  energy  reaching  the  sensor  of  an  electro-opti¬ 
cal  system  or  reduces  the  target/background 
contrast,  or  both.  As  a  result,  many  electro-optical 
systems  become  ineffective  even  during  periods 
of  light  snowfall.  Investigations,  of  the  physical 
and  optical  properties  of  falling  snow  are  neces¬ 
sary  to  evaluate  the  performance  degradation  of 
electro-optical  devices  in  a  winter  environment. 

T  ie  extinction  coefficient  is  a  measure  of  the 
attenuation  of  radiation  as  it  propagates  through 
the  atmosphere.  The  effects  of  falling  snow  on  the 
attenuation  of  optical  radiation  were  investigated 
during  a  series  of  winter  experiments  (SNOW 
tests).  Extensive  transmission  measurements  at 
several  infrared  and  visible  wavelengths  were 
made  in  conjunction  with  detailed  measurements 
characterizing  the  properties  of  falling  snow.  These 
tests  have  led  to  the  refinement  of  empirical  and 
theoretical  models  for  predicting  the  extinction  of 
visible  and  infrared  radiation  by  falling  snow  and 
also  have  led  to  improved  techniques  for  charac¬ 
terizing  the  properties  of  snowfall.  Although  much 
progress  was  made  during  the  SNOW  tests,  our 
knowledge  of  the  physical  and  optical  properties 
of  falling  snow  remains  incomplete,  litis  report 
reviews  what  is  known  and  is  intended  to  provide 
guidelines  for  future  research,  and  also  to  provide 
engineering  and  design  guidelines  for  developing 
efficient  methods  to  predict  and  improve  the  per¬ 
formance  of  electro-optical  systems  in  a  snow 
environment. 


PROPERTIES  OF  FALLING  SNOW 

The  major  obstacle  in  the  study  of  optical  energy 
transmission  through  falling  snow  is  the  complex 
and  diverse  properties  of  snow  particles.  Snow¬ 
storms  occur  under  diverse  meteorological  condi¬ 
tions  so  that  snow  particles  of  different  size  shape 
and  amount  are  generated.  These  snowfall  prop¬ 
erties  must  be  known  in  order  to  develop  models 
for  predicting  the  transmission  of  visible  and  in¬ 
frared  radiation  through  falling  snow.  The  proper¬ 
ties  of  falling  snow  observed  during  the  SNOW 
tests  and  the  techniques  used  to  measure  these 
properties  are  presented. 

Snow  crystal  type 

During  their  initial  growth  phase,  snow  par¬ 
ticles  possess  hexagonal  plane  and  prismatic 
shapes.  These  hexagonal  forms  are  maintained  as 
long  as  the  growth  takes  place  by  vapor  diffusion. 
The  shapes  of  these  crystals  'hat  grow  strictly  by 
diffusion  are  dependent  on  ;he  ambient  tempera¬ 
ture  and  humidity  (Table  1).  As  these  particles 
grow  and  become  large  enough  to  precipitate, 
nondiffusive  growth  mechanisms  such  as  riming 
(accumulation  of  cloud  droplets)  and  aggregation 
produce  more  complicated  snow  structures.  Frag¬ 
mented  snow  particles  are  frequently  observed, 
which  suggests  that  the  snow  particles  are  also 
modified  by  collision  or  other  breakup  mecha¬ 
nisms.  Hie  resulting  variety  of  snow  particle  shapes 
is  evident  from  the  schematic  representation  of 
snow  types  introduced  by  Magono  and  Lee  (1 966), 
shown  in  Figure  1. 

The  detailed  snow  classification  of  Magono  and 
Lee  (1966)  needs  to  be  simplified  for  practical  use. 
At  present,  a  satisfactory  snow  classification 
scheme  useful  for  transmission  studies  has  not 
been  developed.  A  reduced  set  of  snow  types 
should  be  devised  based  on  similar  physical  or 
optical  properties.  A  simple  classification  scheme 
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Figure  2,  Schematic  representation  of  various  forms  of  snow  and  frozen  particles  (after  Magono  and  Lee  1966)., 


Table  1.  Snow  crystal  habit  as  a  function  of  am¬ 
bient  conditions  (°C)  (Houghton  1986). 

At  or  above  Between  ice  and 


water  saturation _ water  saturation 


Oto—4 

plates 

0  to  -4 

plates 

-4  to  -6 

needles 

-4  to -10 

columns 

-6  to  -10 
-9  to -13 
-16  to —40 
-13  to -17 
-17  to -21 
-21  to -40 

columns  -10  to  -20 

plates,  sector  plates 
columns 

stellars,  dendrites 
thick  and  sector  plates 
columns,  plates,  bullets 

thick  plates 

describing  the  degrees  of  riming,  aggregation  and 
deviation  from  the  spherical  shape  may  be  a  useful 
starring  point. 

Tire  Formvar  replication  technique  (Schaefer 
1941)  is  frequently  used  to  analyze  the  shapes  of 
snow  particles.  By  use  of  this  technique,  plastic 
replicas  of  snow  particles  are  made  so  that  snow 
types  can  be  determined  at  a  later  time  in  a  labora¬ 
tory.  A  replica  is  made  by  allowing  a  snow  particle 
to  fall  on  a  surface  that  has  been  previously  coated 
with  a  thin  layer  of  replicating  solution.  The  solu¬ 
tion  is  made  by  dissolving  plastic  in  a  solvent 
(polyvinyl  formal,  better  known  as  Formvar,  dis- 
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Figure  2.  Top  and  front  view  of  a  continuous  snow  replicator  showing  the  control  paneland 
the  film  spools.  Snow  replicas  are  made  on  a  clear  70-mm-wide  motion  film  coated  with 
Fornwar  solution. 


solved  in  1,2-dichloroethane,  is  commonly  used). 
After  a  snow  particle  lands  on  the  coated  surface, 
the  solution  covers  the  snow  through  surface  at¬ 
traction.  The  solvent  then  evaporates,  leaving  the 
snow  particle  encased  in  a  plastic  film.  After  the 
snow  sublimates,  a  hollow  plastic  replica  of  the 
original  crystal  remains. 

Snow  types  can  change,  sometimes  rapidly,  dur¬ 
ing  a  snowfall.  A  continuous  record  of  the  snow 
type  throughout  a  snowstorm  can  be  maintained 


by  automating  the  Formvar  replication  technique. 
Figure  2  shows  a  continuous  replicator  that  was 
used  during  the  SNOW  tests.  A  clear  70-mm-wide 
motion  film  that  has  been  coated  with  a  thin  layer 
of  the  Formvar  solution  moves  and  captures  the 
snow  particles  as  they  fall  through  a  sampling 
opening.  The  film  is  then  wound  onto  a  takeup  reel 
after  it  has  traveled  a  sufficient  distance  to  ensure 
that  the  solvent  has  completely  evaporated.  A  con¬ 
tinuous  length  of  film  with  the  snow  replicas  is 
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produced  that  can  be  analyzed  with  any  suitable 
magnifying  instrument. 

The  analysis  of  snow  types  using  the  Formvar 
replication  technique  is  tedious  and  time-consum¬ 
ing.  Alternative  methods  for  determining  the  snow 
types  are  desirable.  Efforts  are  underway  to  relate 
the  snowfall  properties  at  ground  level  with  spe¬ 
cific  synoptic  meteorological  conditions.  For  exam¬ 
ple,  it  may  be  possible  to  relate  riming  and  aggre¬ 
gation  to  large-scale  weather  patterns.  However, 
at  present,  the  investigations  of  such  relationships 
are  limited  and  the  numerical  models  describing 
precipitation  formation  are  too  simplistic  to  pre¬ 
dict  the  detailed  snowfall  properties  at  specified 
locations. 

Snow  size  distribution 

The  size  distribution  of  snow  particles  is  needed 
to  calculate  the  attenuation  of  optical  energy  by 
falling  snow.  A  widely  used  expression  to  model 
the  snow  size  is  the  Gunn-Marshall  distribution. 
Gunn  and  Marshall  (1958)  segregated  their  snow¬ 
storm  data  according  to  precipitation  rate  and 
used  an  exponential  distribution  to  fit  the  results. 
The  Gunn-Marshall  distribution  is  expressed  as 

n(Dm)  =  nm  exp(— A.m  ■  Dm)  (1) 

where  the  distribution  parameters 

Xm  =  25.5  •  Rf  48  ■  mm'1 

and 

nm  =  3.8  •  103  •  Ks"0  87  nr3  mm'1 

Dm  is  the  diameter  of  a  water  drop  to  which  a  snow 
particle  melts  and  Rs  is  the  precipitation  rate  in 
millimeters  of  water  per  hour.  Sekhon  and  Srivas- 
tava  (1970)  reanalyzed  the  data  of  Gunn  and 
Marshall  and  of  others  and  proposed  that  more 
consistent  results  could  be  achieved  if  the  distribu¬ 
tional  parameters  had  the  values 

Xm  =  22.9  ■  Rs-0  54  mm"1 


and 

«m  =  2.5  •  103  •  R^m  mr3  mm'1 . 

The  simplicity  of  the  exponential  distribution 
makes  it  ideal  for  use  in  theoretical  studies.  The 
moments  i  of  the  exponential  distribution  are  given 
by 


1  \  #  r-»(l+l) 

Dm  M(  Dm)  dD  m  =  M m  1  ~ 


where  T  is  the  gamma  function.  This  equation 
implies  that  there  is  a  correlation  between  the 
average  snow  size  (first  moment  of  the  distribu¬ 
tion)  and  the  snow  precipitation  rate.  However, 
the  dependence  of  snow  size  distribution  on  pre¬ 
cipitation  rate  alone  may  be  an  oversimplification 
when  one  considers  the  large  variability  in  the 
structure  of  the  snow  particles.  The  observations 
of  Feng  and  Grant  (1982)  on  the  correlation  of  the 
snow  type,  the  number  flux  and  the  snow  intensity 
indicate  that  the  snow  size  is  not  necessarily  de¬ 
pendent  on  the  precipitation  rate.  The  results  ob¬ 
tained  during  the  SNOW  tests  also  indicate  that 
the  snow  size  is  independent  of  the  snow  precipi¬ 
tation  rate. 

During  the  SNOW  tests,  100  to  200  snow  par¬ 
ticles  were  replicated  on  10  by  13  cm  glass  plates  at 
the  start  of  each  hour  during  a  snowstorm.  The 
areas  of  the  snow  particles  were  measured  by 
projecting  the  video  images  of  the  snow  particles 
on  a  TV  monitor  and  outlining  them  with  a  light 
pen.  A  computer  was  used  to  measure  the  outlined 
area  of  each  particle  and  to  compute  its  equivalent 
radius  (the  radius  of  a  circle  with  the  same  area). 
Using  this  method  the  average  snow  sizes  were 
obtained  and  compared  with  the  snow  mass  con¬ 
centration.  Figure  3  shows  that  a  correlation  be¬ 
tween  snow  size  and  snow  mass  concentration 
was  not  observed  during  the  SNOW  tests,  which 
contradicts  the  exponential  size  distribution  model 


Figure  3.  Relationship  between  the  average  snow  size  and 
the  snow  mass  concentration.  A  correlation  between  snow 
size  and  snow  intensity  does  not  exist. 
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(although  snow  mass  concentration  is  used  in¬ 
stead  of  snow  precipitation  rate,  the  argument  re¬ 
mains  valid  since  a  correlation  exists  between  the 
two  snowfall  parameters).  At  the  lower  snowfall 
intensities,  the  average  snow  size  ranged  between 
0.3  and  1 . 1  mm  in  radius.  However,  as  the  snowfall 
intensity  increased,  the  average  snow  size  was  less 
than  0.8  mm. 

The  snow  size  distributions  observed  during 
the  SNOW  tests  varied  significantly  and  it  was  dif¬ 
ficult  to  find  suitable  analytical  expressions  to 
model  the  snow  size.  The  collected  snow  samples 
frequently  contained  several  snow  types  with  vary¬ 
ing  degrees  of  riming  and  aggregation,  so  the  data 
could  not  be  easily  segregated  according  to  snow 
type.  Therefore,  the  snowstorms  were  segregated 
based  on  the  average  snow  size,  and  the  modified 
gamma  distribution  was  used  to  describe  the  snow 
size  distribution  for  each  class  of  snowstorm.  The 
modified  gamma  distribution  is  a  three-parameter 
model  expressed  as 

n(a )  =  k  ■  a  “•  expK-a/yXrt  /rtc)Y]  (3) 

where  n(a)  is  a  continuous  and  an  integrable  func¬ 
tion  of  radius  a.  The  shape  of  the  distribution  is 
governed  by  a  (positive  integer),  y  (positive  con¬ 
stant)  and  ac  (mode  radius).  The  value  K  is  used  to 
normalize  the  total  number  of  particles  in  a  unit 
volume  and  does  not  affect  the  shape  of  the  distri¬ 
bution. 

Figure  4  illustrates  the  modified  gamma  distri¬ 
bution  fit  for  the  composite  size  distributions  for 


snowfalls  whose  average  radii  range  from  0.2-0.7 
mm.  These  size  ranges  accounted  for  approxi¬ 
mately  90%  of  the  snowfalls  observed  during  the 
SNOW  tests.  These  size  distributions  were  ob¬ 
tained  from  snow  particles  that  were  replicated  on 
a  horizontal  surface.  To  infer  the  size  distribution 
in  space  (unit  volume),  the  observed  distributions 
must  be  divided  by  the  fall  velocity  of  the  snow 
particles  in  each  size  interval.  Since  the  snow  fall 
velocity  is  not  uniquely  dependent  on  the  size, 
errors  are  introduced  when  inferring  the  size  dis¬ 
tribution  in  space  from  the  replication  technique. 

The  segregation  of  the  snow  data  based  on  the 
average  size  only  provides  a  rough  estimate  of  the 
actual  distribution.  The  composite  distributions 
do  not  adequately  represent  the  snow  size  distri¬ 
bution  for  snowfall  that  is  composed  of  large 
aggregates.  A  bimodal  distribution  is  often  ob¬ 
served  as  the  number  of  the  large  aggregates 
increases.  To  more  accurately  represent  the  snow 
size  distributions,  an  analytical  expression  that 
takes  into  account  the  bimodal  distribution  is 
required.  At  present,  information  on  the  size  dis¬ 
tribution  of  snow  particles  is  limited  and  more 
studies  are  required  to  develop  realistic  models  of 
the  snow  size. 

To  obtain  more  information  about  the  snow  size 
distribution,  improved  techniques  for  measuring 
snow  particle  sizes  are  needed.  The  available  in¬ 
formation  on  the  snow  sizes  has,  for  the  most  part, 
been  obtained  using  laborious  methods  such  as 
replication  or  photography.  To  improve  the  snow 
measurement  capability,  Berger  (1983)  investi- 


a.  Average  size  02-0.3  mm  (a  =  5,y=  0.5,  b.  Average  size  0.3-0.4  mm  (a~4,y=  0.S, 

a  =  0.15),  a  =  0.25). 

*  t 

Figure  4.  Modified  gamma  distribution  fit  for  composite  snow  size  distributions.  The  data  went  segre¬ 
gated  according  to  the  average  snow  size.. 
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c.  Average  size  0.4-0. 5  mm  (a  =  3,  y=  1.1, 
ac  =  0.32). 


e.  Average  size  0.6-0.7  mm  (a=  1,  y=  1.9, 
ac  =  0.45). 


Figure  4  (cant'd).  Modified  gamma  distribu¬ 
tion  fit  for  composite  snow  size  distributions. 
The  data  were  segregated  according  to  the  aver¬ 
age  snow  size. 

gated  the  feasibility  of  vising  a  PMS  (Particle 
Measuring  System)  two-dimensional  optical  ar¬ 
ray  imaging  probe  to  measure  the  snow  sizes  dur¬ 
ing  the  SNOW  tests.  The  imaging  probe  uses  a 
linear  array  of  photodiodes  as  a  size-measuring 
grid  in  a  shadowgraph-type  manner.  A  particle 
entering  the  sample  volume  passes  through  a  beam 
of  light  that  produces  a  shadow  on  the  photodi¬ 
odes.  The  photodiodes  are  scanned  at  a  known  fre¬ 
quency  to  produce  a  digitized  image  of  each  par¬ 
ticle.  The  data  are  in  a  form  that  can  be  easily  man¬ 
aged  with  computer  processing  so  that  large 
amounts  of  data  can  be  obtained  with  a  minimum 
amount  of  attention.  A  schematic  of  the  probe  and 


d.  Average  size  0.5-0. 6  mm  (a  =  2,  y=  1.0, 
ac  =  0.35). 


some  examples  of  digitized  images  of  snow  par¬ 
ticles  are  shown  in  Figure  5. 

The  PMS  probes  were  designed  to  be  aircraft- 
mounted  to  measure  the  size  of  cloud  and  precipi¬ 
tating  particles.  For  ground  operations,  the  snow 
particles  were  aspirated  into  the  sample  volume  at 
a  known  velocity  (since  the  photodiodes  are 
scanned  at  fixed  frequency,  the  velocity  of  the  par¬ 
ticle  as  it  passes  the  beam  must  be  known).  To  in¬ 
terpret  the  data  obtained  from  ground  operations, 
the  effect  of  the  aspiration  on  the  measurements  of 
number  concentration  and  size  distribution  of  snow 
particles  was  investigated.  The  numerical  studies 
by  Norment  (1986)  and  the  experimental  studies 
by  Humphries  (1985)  and  Holroyd  (1986)  on  the 
effects  of  aspiration  are  not  conclusive  and  more 
studies  are  required.  When  the  PMS  probes  can  be 
properly  calibrated  for  the  effect  of  aspiration, 
their  ability  to  measure  the  snow  sizes  rapidly  and 
continuously  makes  them  ideal  instruments  for 
future  snow  size  distribution  studies. 

Snow  size  versus  ambient  conditions 

The  average  snow  size  (equivalent  radius)  ob¬ 
served  during  the  SNOW  tests  generally  ranged 
from  0.2  to  0.7  mm.  The  frequency  distribution  of 
the  average  size  is  shown  in  Figure  6.  It  was 
previously  shown  that  the  average  snow  size  was 
not  related  to  the  snowfall  intensity  (Fig.  3).  The 
relationships  between  the  average  snow  size  and 
the  surface  meteorological  conditions  measured 
during  the  SNOW  tests  were  also  investigated. 
Figures  7  and  8  illustrate  the  relationships  of  the 
average  snow  size  to  surface  temperature  (meas¬ 
ured  at  2-m  height)  and  to  relative  humidity  re¬ 
spectively.  For  snowfalls  that  occurred  when  the 
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Relative  Frequency  (%) 


Beam  Cross  Section 


Figure  5.  Schematic  of  a  Particle  Measuring  System  optical  imaging  probe  and 
digitized  images  of  snow  particles  (from  Berger  1983). 


Average  Snow  Radius  (mm)  Temperature  (°C) 


Figure  6.  Distribution  of  the  average  snow  Figure  7.  Relationship  between  the  average  snow  size  and 

size  (equivalent  radius)  obtained  from  Form -  the  surface  temperature, 

var  replicas. 
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Average  Snow  Radius  (mm) 


1954)  and  arrived  at  the  following  relationship 


r - 1 - 1 - t - 1 - ; 


60  70  80  90  100 

Relative  Himiaity  (%) 


Figure  8.  Relationship  between  the  average  snow  si ze  and 
the  relative  humidity . 


surface  temperature,  T,  was  between  -18and  O'C, 
a  linear  regression  of  the  data  resulted  in  the 
following 

a  (mm)  =  0.6  +  0.021 .  v4) 

The  relationship  between  snow  size  and  relative 
humidity  (RH  is  expressed  using  decimals)  re¬ 
sulted  in  the  following  linear  regression 

a  (mm)  =  -0.33+RH.  (5) 

Although  the  figures  indicate  that  the  snow  sizes 
increase  as  the  temperature  and  the  humioity  in¬ 
crease,  the  scatter  in  the  data  is  such  that  the  re¬ 
gression  results  may  not  be  very  meaningful. 

Snow  particle  fall  velocity 
The  fall  velocities  of  snow  particles  are  impor¬ 
tant  when  investigating  the  relationship  between 
snow  precipitation  rate  and  snow  extinction  coef¬ 
ficient.  Due  to  the  diverse  sizes  and  shapes  of  snow 
particles,  a  range  of  snow  fall  velocity  exists.  Nu¬ 
merous  empirical  relationships  between  the  snow 
fall  velocity  and  the  snow  type  and  size  have  been 
reported.  Typically,  the  observed  data  are  fitted  by 
a  power  function  of  the  type 

vs -a -dm  (6) 

where  vs  is  the  fall  velocity,  dw  is  the  size  of  the 
snow  particle  and  a  and  b  are  contants  for  a  given 
type  of  snowfall  (Jiusto  and  Bosworth  1971,  Lo- 
catelli  and  Hobbs  1974). 

Other  empirical  fits  to  the  snow  fall  velocity  data 
have  also  been  reported.  O'Brien  (1970)  analyzed 
thedataof  Nakaya  (1954)  andMagono(1951, 1953, 


Vs  =  a  -dm- Dogin(ys-d m)+b]  (7) 

where  a,  b,  and  r.  are  empirical  constants  that  were 
chosen  to  fit  the  data.  Suzuki  et  al.  (1981)  studied 
the  relationship  between  the  snow  fall  velocity 
and  snow  density,  ps(g/cm3),  and  observed  the 
following: 


‘2.36-  p P-4 


(8) 


b.  vs  =2.36  ■  pF  •  dm-23 (after  Suzuki  et  al.  1981). 


Figure  9.  Empirical  relationships  between  snowfall 
velocity  and  size  for  different  snow  types. 
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Average  Snow  Fall  Velocity  (m/s) 

b.  Rosinski  et  al.  (1983).. 


Figure  10.  Distribution  of  the  average  snoiu  fail  velocities. 


Since  the  methods  used  to  describe  the  size  and  to 
classify  the  snow  type  often  vary  among  research¬ 
ers,  it  is  difficult  to  compare  the  results  obtained  by 
the  different  investigators,  figure  9  illustrates  the 
empirical  fits  for  different  types  of  snowfall  re¬ 
ported  by  O'Bnen  (1970)  and  Suzuki  et  al.  (1981). 

Meilor  (1 966)  and  Rosinski  et  al.  (1983)  observed 
the  fall  velocities  of  snow  particles  over  a  period 
of  few  minutes  to  find  the  average  snow  fall  veloc¬ 
ity.  The  average  fall  velocity  reported  by  these 
investigators  generally  ranged  from  0.4  to  1.1  m/ 
s.  The  frequency  distributions  of  the  average  fall 
velocities  are  shown  in  Figure  10.  The  results  are 
very  similar,  and  both  distributions  have  the  mode 
value  of  0.75  m/s  for  the  fall  velocity. 

During  the  SNOW  tests,  Berthel  et  al.  (1983) 
measured  snow  particle  fall  velocities.  Multiple 
images  of  a  snow  particle  were  recorded  on  a  sin¬ 


gle  video  frame  by  illuminating  the  vertical  fall 
with  a  strobe  light.  The  fall  velocity  of  individual 
particles  was  determined  from  the  distance  be¬ 
tween  the  images  and  the  known  frequency  of  the 
strobe  illumination.  A  power  fit  was  used  to  relate 
the  snow  size  and  the  fall  velocity.  They  found  that 
the  snow  fall  velocities  generally  ranged  between 
0.5  and  2.0  m/s,  which  agreed  well  with  the  previ¬ 
ously  described  empirical  relationships. 

Snow  particle  mass 

The  mass  (or  density)  of  individual  snow  parti¬ 
cles  is  difficult  to  measure  and  large  discrepancies 
exist  among  various  investigators.  The  relation¬ 
ships  of  snow  precipitation  and  snow  mass  con¬ 
centration  to  the  optical  attenuation  is  greatly 
affected  by  the  mass  of  the  snow  particles.  Figure 
11  shows  the  range  of  snow  masses  observed  by 
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Figure  11.  Range  of  snow  masses  for  commonly  observed 
snow  types  (after  Locatelli  and  Hobbs  1974). 
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Average  Snow  Mass  (mg) 

Figure  12.  Distribution  of  the  average  snow 
mass  (after  Rosinski  et  al.  1983). 

Locatelli  and  Hobbs  (1974).  The  snow  mass  was 
determined  by  melting  the  snow  particles  in  oil 
and  measuring  the  melted  droplet  diameter. 
Rosinski  et  al.  (1983)  measured  the  average  mass 
of  snow  particles  during  snowfall  by  both  con¬ 
ductive  titration  and  by  weighing  the  individual 
particles.  A  histogram  of  the  average  mass  ob¬ 
served  by  these  investigators  is  shown  in  Figure 
12.  No  attempt  was  made  to  measure  the  mass  of 
individual  snow  particles  during  the  SNOW  “ests. 

Snow  precipitation  rate 
The  most  common  quantitative  measure  of  the 
falling  snow  is  the  precipitation  rate.  This  is  typi¬ 
cally  measured  by  monitoring  the  amount  of  snow 
that  falls  through  a  horizontal  surface  sample  area 


mto  a  tippmg-Ducket  or  weighing-type  device. 
These  mechanical  gauges  have  poor  precision  and 
time  resolution,  so  they  are  inadequate  to  support 
electro-optical  system  performance  tests  where 
the  short-term  variations  in  the  snowfall  rate  are 
important.  Berthel  et  al.  (1983)  and  Lacombe  (1983) 
used  electronic  balances  to  improve  the  weighing 
technique  for  measuring  the  snow  precipitation 
rate.  They  were  able  to  reduce  the  average  sam¬ 
pling  time  for  measuring  precipitation  rate  to  1 
minute. 

An  optical  technique  for  measuring  snowfall 
rate,  which  may  be  more  suitable  than  the  elec¬ 
tronic  balance  gauge  foi  field  applications,  was  in¬ 
vestigated  by  Koh  and  Lacombe  (1986).  The  opti¬ 
cal  gauge  shown  in  Figure  13,  which  was  origi¬ 
nally  designed  to  measure  rain  rates  (Wang  et  al. 
1980),  was  modified  for  operation  in  snow.  The 
light  source  for  the  optical  gauge  is  an  infrared- 
emitting  diode  that  is  partially  collimated.  Snow 
particles  falling  through  the  sample  area  break  the 
light  beam  to  produce  light  intensity  fluctuations 
on  the  receiver.  The  receiver  output  is  passed 
through  a  bandpass  filter  (75-900  Hz),  and  its  rms 
value  is  measured.  This  yields  a  slowly  varying 
(10-s  time  constant)  dc  voltage  proportional  to  the 
snow  rate.  Comparisons  of  the  precipitation  rates 
measured  with  an  electronic  balance  gauge  an  d  an 
optical  gauge  have  shown  that  the  snc.w  precipita¬ 
tion  rate  ran  be  estimated  within  a  factor  of  two 
using  the  optical  gauge.  The  loss  of  accuracy  in 
snow  rate  measurement  using  the  optical  tech¬ 
nique  is  compensated  foi  by  its  simplicity  (virtu¬ 
ally  no  maintenance  is  required)  and  improved 


Figure  13.  Optical  gauge  for  measuring  snowpredpitation  rate. 
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time  resolution  over  conventional  snow  gauges. 
The  accuracy  of  snow  rate  measurement  can  be 
improved  by  segregating  the  data  according  to 
snow  types.  Studies  are  being  conducted  to  deter¬ 
mine  the  feasibility  of  identifying  the  snow  types 
using  the  optical  gauge. 

Snow  mass  concentration 

The  snow  mass  concentration  is  the  mass  of 
snow  per  unit  volume  of  the  atmosphere  during  a 
snowstorm.  Stallabrass  (1976)  introduced  a  direct 
means  of  measuring  the  snow  mass  concentration 
that  was  later  improved  by  Lacombe  (1983).  Fig¬ 
ure  14  shows  the  device  used  to  measure  the  snow 
mass  concentration  during  the  SNOW  tests.  A  col¬ 
lection  head  mounted  on  an  arm  rotating  in  a 
horizontal  plane  sweeps  a  known  volume  of  air 
and  the  air-snow  mixture  enters  a  collection  ori¬ 
fice.  The  snow  particles  are  then  separated  from 
the  air,  melted  and  forced  through  a  needle  by 
centrifugal  force.  Equally  sized  droplets  form  at 
the  needle  tip  and  are  photoelectricaliy  counted  as 
they  separate  from  the  needle  and  pass  between  an 
infrared-emitting  diode  and  a  phototransistor.  The 
mass  of  the  droplet  is  calculated  from  a  force- 
balance  analysis  that  relates  centrifugal  and  sur¬ 
face  tension  forces.  Summing  the  nu  mber  of  drop¬ 
lets  over  a  given  time  period  provides  a  measure  of 
the  snow  mass  concentration. 

Accurate  measurement  of  snow  mass  concen¬ 
tration  requires  a  high  collection  efficiency  rate  for 
the  collection  inlet.  A  numerical  study  by  Norment 
(1986)  determined  the  sampling  efficiency  of  the 


mass  concentration  measuring  instruments  to  be 
excellent  over  a  compreheiisive  size  range  of  water 
drops  (20-  to  1000-pm  diameter)  and  hexagonal  ice 
plates  (40-  to  700  -pm  plate  diameter).  The  opera¬ 
tion  of  these  instruments  during  the  SNOW  tests 
demonstrated  that  the  snow  mass  concentration 
can  be  measured  over  extended  time  periods  and 
that  the  changes  in  snow  intensity  can  be  rapidly 
monitored  (once  every  minute),  making  this  in¬ 
strument  a  valuable  tool  for  characterizing  the 
snow  environment  during  tests  of  electro-optical 
system  performance. 

Relationship  between  mass  concentration 
and  precipitation  rate 

The  relationship  between  the  snow  mass 
concentration  Ms  and  the  snow  precipitation  rate 
Rs  can  be  approximated  by 


Ms  —  Rs  /  Vs 


(9) 


where  Fis  some  representative  snow  fall  velocity. 
Since  Rs  is  commonly  measured,  the  ability  to 
estimate  Ms  from  Rs  measurement  is  of  practical 
interest.  Mellor  (1966, 1983)  measured  Rsand  the 
average  snow  fall  velocity  to  estimate  Ms.  The 
relationship  between  Ms( g/m3)  and  Rs  (mm/hr) 
was  described  by  the  following  relationship 


Ms  =  0.32-  Rs 


0.9 


(10) 


Koh  et  al.  (1988)  directly  measured  Msand  Rsand 
the  linear  regression  fit  to  the  data  resulted  in  the 
following  Ms  to  Rs  relationships 


40 


Figure  15.  Relationshipbetween  the  snow  precipita¬ 
tion  rate  and  the  snow  mass  concentration. 


Ms  =  0.30  Rs  ("dense"  snow)  (11a) 

Ms  =  0.47- Rs  ("light"  snow) .  (lib) 

The  "dense"  and  "light"  classifications  refer  to  the 
bulk  density  of  snow  crystals.  Snowstorms  com¬ 
posed  of  ice  pellets,  needles,  columns  and  rimed 
particles  were  classified  as  "dense"  snow.  Snow 
particles  with  cavities  and  irregularities  in  the 
external  crystal  structure,  such  as  aggregates  of 
snow  particles  and  spatial  dendrites,  were  classi¬ 
fied  as  "light"  snow.  Although  the  classification  is 
subjective,  the  Ms  to  Rs  relationships  for  the  two 
types  differ.  Figure  15  illustrates  the  Ms  to  Rs 
relationships  observed  by  Mellor  (1983)  and  Koh 
et  al.  (1988). 


THEORY  OF  EXTINCTION 


The  amount  of  transmitted  optical  energy  re¬ 
ceived  by  a  detector  located  downstream  is  re¬ 
duced  during  a  snowstorm.  This  process,  called 
extinction,  is  caused  by  the  energy  redistribution 
of  the  transmitted  radiation  due  to  scattering  and 
absorption  by  the  snow  particles  along  the  trans¬ 
mission  path.  The  extinction  characteristic  of  snow 
particles  is  described  in  terms  of  the  extinction 
coefficient  Bwt,  which  is  expressed  as 


Rex  t  — 


[Csca(fl  )+Cabs(a )]  n{a)da.  (12a) 


Jo 


C^and  Cabs  represent  the  scattering  and  the  ab¬ 
sorption  cross  sections,  respectively,  and  «(a)  is 
the  number  of  particles  per  unit  volume  with  sizes 
between  a  and  a+da.  The  total  energy  scattered  in 


all  directions  and  the  total  energy  absorbed  inside 
a  particle  are  equal  to  the  energy  of  the  incident 
wave  falling  on  the  areas  C^and  (^respectively. 
The  extinction  cross  section  C  is  defined  as 

cs„+cabs. 

The  dimensions  of  the  snow  particles  are  much 
larger  than  the  visible  and  infrared  wavelengths; 
therefore,  short-wavelength  limit  approximation 
can  be  used  to  determine  the  extinction  cross  sec¬ 
tion.  In  this  approximation,  the  extinction  cross 
section  is  equal  to  twice  the  geometrical  cross 
section  (van  de  Hulst  1957).  Equation  12a  can 
therefore  be  simplified  to 


Bext 


jc  a 2  n(a)  da 


Jo 


(12b) 


where  it a}  is  the  area  of  a  snow  particle  projected 
normal  to  the  incident  beam.  This  seeming  contra¬ 
diction,  that  a  large  particle  removes  from  the 
incident  radiation  twice  the  amount  of  energy  it 
can  intercept,  is  referred  to  as  the  extinction  para¬ 
dox. 

The  extinction  paradox  can  be  explained  using 
geometrical  optics  and  physical  optics  approx¬ 
imations.  The  geometrical  optics  approximation 
(subdividing  the  incident  radiation  into  a  large 
number  of  rays  and  using  the  Fresnel  equation  and 
Snell's  law  to  determine  the  interaction  at  the 
boundaries)  indicates  that  the  energy  incident  on 
the  geometrical  cross  section  m2  is  either  reflected 
or  refracted.  The  refracted  energy  is  either  ab¬ 
sorbed  or  it  is  refracted  out  of  the  particle.  In 
addition  to  this  reflected  and  refracted  cross  section 
the  diffraction  also  contributes  to  the  total  ex¬ 
tinction. 

The  effects  of  diffraction  can  not  be  explained 
using  geometrical  optics;  therefore,  a  physical  op¬ 
tics  solution  (scalar  diffraction  theory)  is  required. 
Invoking  the  "Babinet  Principle,"  which  states 
that  the  diffraction  patterns  are  similar  when  ob¬ 
served  behind  either  an  aperture  or  an  obstacle  of 
the  same  shape  and  size,  one  can  show  that  the 
amount  of  diffracted  energy  is  also  equal  to  the 
amount  of  energy  incident  on  the  geometrical 
cross  section  m2.  The  combined  effects  of  reflec¬ 
tion,  refraction  and  diffraction  result  in  the  extinc¬ 
tion  cross  section  of  2m2. 

Once  the  extinction  coefficient  is  known,  the 
variation  cf  the  radiation  intensity  1  along  the 
propagation  path  L  can  be  determined  using  the 
relationship 


dl  =  —/  Bext  ■  dL  . 


(13) 
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Figure  16.  Illustration  of  the  forward-scattered 
light  measured  by  a  detector.  Measured  extinction 
is  less  than  the  theoretical  by  the  amount  of  the 
forward-scattered  light  contained  in  the  shaded 
region.. 


Integration  of  this  equation  yields  the  familiar 
Beer-Lambert  relationship  for  the  "theoretical" 
transmittance  T  (the  reason  for  specifying  "theo¬ 
retical"  is  explained  later)  through  the  atmosphere 

T  =  i  =  exp(-Be«-  L)  (14) 

Jo 

where  I„is  the  intensity  at  the  start  of  the  path  and 
I  is  the  intensity  at  distance  L. 

Equation  14  is  valid  if  the  contribution  to  the 
received  power  from  forward  scattering  is  negli¬ 
gible.  In  actual  practice,  since  all  detectors  have  a 
nonzero  field  of  view,  a  portion  of  the  radiation 
scattered  in  the  forward  direction  contributes  to 
the  received  power.  The  effect  of  forward  scatter¬ 
ing  is  illustrated  in  Figure  16.  The  shaded  area 
represents  the  forward-scattered  radiation  enter¬ 
ing  the  acceptance  angle  of  a  detector.  The 
"theoretical"  prediction  for  transmittance  is  al¬ 
ways  less  than  the  "measured"  transmittance,  since 
eq  14  does  not  account  for  forward  scattering 
contribution  to  the  measured  transmittance.  The 
difference  between  the  "measured"  and  the  "theo¬ 
retical"  values  is  determined  by  the  detector  size 
and  the  amount  of  forward-scattered  radiation. 

This  brief  overview  has  been  presented  to  show 
that  the  extinction  coefficient,  the  angular  distri¬ 
bution  of  scattered  radiation,  and  the  geometries 
of  the  transmitted  beam  and  the  receiver  are  re¬ 
quired  to  predict  the  amount  of  radiation  meas¬ 
ured  by  any  electro-optical  device. 

Comparison  of  visible  and 
infrared  transmission 

The  transmission  data  obtained  during  the 
SNOW  tests  indicate  that  snow  is  slightly  more 
transparent  to  visible  light  than  to  infrared  radia¬ 
tion.  Typical  visible  and  infrared  transmission 
data  through  the  falling  snow  are  shown  in  Figure 
17.  Transmission  in  the  far-infrared  (8-  to  14-pm) 
and  mid-infrared  (3-  to  5-pm)  regions  is  lower  than 
in  the  visible  regions.  The  transmission  at  the  near- 
infrared  (1 .06  pm)  region  is  slightly  higher  than  at 
the  visible  wavelength  (in  most  cases,  visible  trans¬ 
mission  through  snowfall  is  better  than  1.06  pm 
regions).  These  measured  transmission  results  are 


contrary  to  eq  12b,  which  indicates  that  extinction 
by  snow  particles  is  independent  of  wavelength  in 
the  visible  and  infrared  regions.  Table  2  shows 
some  of  the  previously  reported  relationships 
between  infrared  and  visible  transmission. 

The  observed  wavelength  dependence  appears 
to  be  the  result  of  the  increase  in  the  scattering  in 
the  near  forward  direction  at  the  shorter  wave¬ 
lengths.  As  the  ratio  of  the  particle  size  to  the 
wavelength  increases,  more  of  the  scattered  radia¬ 
tion  is  concentrated  in  the  forward  direction. 
Therefore,  the  contribution  of  the  narrow-angle 
foiward  scattering  into  the  receiver  is  greater  at 
the  shorter  wavelengths.  The  effect  of  forward 
scattering  on  the  measured  transmission  is  dis¬ 
cussed  in  the  next  section. 

Forward-scattering  c  orrection 

The  angular  distribution  of  the  light  scattered  by 
irregularly  shaped  particles  is  not  known.  There¬ 
fore,  calculations  of  the  forward-scattering  correc¬ 
tion  to  extinction  in  falling  snow  have  for  the  most 
part  been  based  on  the  theories  of  scattering  by 
spherical  particles.  Miller  (1978),  Mill  and  Shettle 
(1983)  and  Seagraves  (1983)  used  spherical  ap-. 
proximation  for  snow  particles  and  used  the  Fraun¬ 
hofer  diffraction  theory  to  calculate  the  intensity 


Time  of  Day 

Figu  re  1 7.  Comparison  of  the  visible  and  infrared 
transmission  data  through  falling  snow.  Snow  is 
slightly  more  transparent  to  visible  light  than  to 
infrared  radiation. 
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Table  2.  Relationships  between  visible  and  infrared  extinction  in  tailing 
snow  (after  Seagraves  1986). 


8"xt(10  6  Jim)  —  1.4-6ext(063  nm* 


£v\t(  lll.b)  —1.38-  Bext(0.63) 

Bext(3  0-5.0)  =1.05-  bJx?(0?5) 

0  09 

Bext(8  0-12  0)  =!.3C  Bext(0.55) 

0  9(1 

Bext(8.0-12.U)  ”  '  25  Be.ir<.  isible) 
B?xt(104)=l-42-  f>ext(0.55) 

Bexttf  0)  =1 .43  •  Bext(0.55/ 

^ext(1.06)  =1-06'  ®ext(0.55) 
Bext(l 0.5'  =1-18'  8ext(0.55) 


of  light  scattered  in  the  near-forward  direction. 
Using  the  Fraunhoher  diffraction  theory,  the  in¬ 
tensity  in  the  forward  direction  can  be  expressed 
as 


dCsea  _  X4  Ji(ssine) 
dCl  k2  xsinG  . 


(15) 


where  the  differential  cross  section  dC dQ.  is  the 
amount  of  radiation  scattered  per  unit  incident  ra¬ 
diation  into  a  unit  solid  angle  about  a  given 
direction,  k  is  the  wave  number  (2k/ X),  x  is  the  size 
parameter  (fai),  and  /,  is  the  first-order  Bessel 
function.  0  is  the  scattering  angle,  and  0  and  180* 
represent  the  forward  and  backward  scattering 
angles  respectively.  A  normalized  scattering  dia¬ 
gram  as  a  function  of  (x  sinG)  is  shown  in  Figure 
18. 

The  amount  of  the  forward-scattered  ligh  t  enter¬ 
ing  a  detector's  field  of  view  of  can  be  calculated  by 
integrating  eq  15  over  the  acceptance  angle  of  the 
detector.  For  a  circular  detector  i  hat  subtends  half¬ 
angle  0  at  a  particle  =  29),  the  equation  be¬ 
comes 


2k  x4  l  j  M*  sine)  l2 

k2  I  L  xsinG  J 


sinG  d  0. 


(16) 


Integrating  eq  16  and  then  normalizing,  one  ob¬ 
tains  the  following 


F  =[l -0  (x  sinG)  -  /?  ( x  sinG)]  (17) 

where  F  is  the  relative  amount  of  the  diffracted 
light  collected  by  a  detector,  and  /„and  /,are  Bessel 


Bisyarin  etal.  (1971) 

Bisyarinetal.  0973 

Duncan  (1981);  Sola  and  Bergmann  (1977) 


lacombe  and  O  Brien  (1982) 


Seagraves  and  Ebersole  (1983) 


Figure  IS.  Normalized  scattering  in  the  near- for¬ 
ward  direction  for  a  spherical  particle  calculated 
using  the  Fraunhofer  diffraction  theory. 


Figure  19.  Relative  amount  of  the  diffracted  light 
collected  by  a  detector  calculated  as  a  function  of  the 
size  parameter  and  the  angle  subtended  by  a  detec¬ 
tor.  Unear  approximation  is  also  shown. 
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functions  of  order  0  and  1  respectively.  The  value  and  for  oblates, 
of  F  as  a  function  of  (x-sin0)  is  shown  in  Figure  i9. 

w  Qf  —  n 

The  figure  also  illustrates  that  eq  17  can  be  simpli¬ 
fied  by  a  linear  approximation  where  P = (fl  2  sin2©  +  c 2  cos2©)1 11 . 


F  =  0.3  (x  sin6)  xsin0<2.7  (18a) 

F  =  0.8+0.15  (x  sin0)  x  sin0  )  2.7  (18b) 

Since  the  acceptance  angle  of  a  typical  transmis¬ 
sion  system  is  very  small,  the  value  of  F  can,  in 
most  cases,  be  approximated  by  eq  18a. 

The  value  of  F  expressed  in  eq  18  is  appropriate 
for  a  single  particle.  The  calculation  forthe  mean 
value  of  tl:e  fraction  of  diffracted  light  F  collected 
by  a  detector  during  transmission  measurements 
requires  lengthy  integration  since  F  depends  on 
the  size  of  the  snow  particle  and  on  the  location  of 
the  particle  along  the  transmission  path/To  avoid 
the  lengthy  integration,  the  value  for  F  was  ap¬ 
proximated  to  be  equivalent  to  the  value  of  F 
calculated  for  the  mean  particle  size  located  at  the 
midpoint  of  the  transmission  path. 

Since  one-half  of  the  extinction  cross  section  is 
attributable  to  the  diffraction,  the  relationship  be¬ 
tween  the  "measured"  and  the  "theoretical"  ex¬ 
tinction  coefficient  is  given  by 

Bex t  =Bext  (1-0.5  F) 


tb  is  the  angle  that  the  major  axis  makes  with  the 
normal  to  the  plane  onto  which  the  spheroid  is 
projected.  The  size  parameter  x  is k  V a/p  . 

If  the  spheroid  is  randomly  oriented,  the  scatter¬ 
ing  must  be  averaged  over  all  possible  rotations  of 
Lhe  spheroid.  The  average  differential  scattering 
cross  section  then  becomes 

/rfCsca.j  _  f  f  p(Qy<$) -JGcz  s\nQ  d®  (21) 
l da  /  Jo  Jo  dQ. 


where  P(0, 0)  is  the  probability  distribution  of  the 
spheroid  orientation.  These  calculations  indicate 
that  the  fraction  of  the  forward-scattered  light  is 
less  for  nonspherical  particles  than  it  is  for  the 
equivalent  sphere  approximation. 


Precipitation  rate  and  extinction 
The  snow  precipitation  rate  can  be  approxi¬ 
mated  by 


)«(u)  v(a)  m(a)  da 
o  ^ 


(22) 


where  Bex t  is  the  measured  extinction  coefficient. 
As  expected,  B„, approaches  the  theoretical  as  the 
amount  of  the  forward-scattered  light  that  is  col¬ 
lected  by  a  detector  becomes  increasingly  smaller. 

Bohren  and  Koh  (1985)  investigated  the  for- 
ward-scattered  correction  required  for  nonspheri¬ 
cal  particles.  Prolate  and  oblate  spheroids  were 
used  to  model  the  shape  of  the  snow  particles.  A 
spheroid  is  formed  by  rotating  an  ellipse  about  one 
of  its  axes;  therefore,  the  projection  of  a  spheroid 
onto  a  plane  is  an  ellipse.  The  differential  scatter¬ 
ing  cross  section  for  a  spheroid  in  the  near  forward 
direction  is  given  by 


dC^ss.  -x4  fyi  (*7  sing)] 2 
dSl  k 2  (*Y  sin0)  j 

where  y  =[(a/p)  cos2*})  +  ({Vaj  sin2<j>]1/2. 

If  a  and  c  are  the  semimajor  and  the  semimi¬ 
nor  axes  of  the  spheroid,  the  semimajor  a  and 
the  semiminor  p  axis  lengths  of  the  correspond¬ 
ing  projected  ellipses  for  prolates  are 

a  =(a  2sin20  +c2cos?0)1/2 
P  =  c 


where  v(a)  =  fall  velocity 
m(a)  =  mass  of  snow 
pw  =  density  of  water. 


Combining  eq  12  and  22,  one  can  express  the  rela¬ 
tionship  between  theextinction  coefficient  and  the 
snow  precipitation  rate  as 

B?xt=k  •  R$  (23) 

where  k  is  equal  to 


271  a  2  ti{a )  da 


k  -  pw — - 


n(a)  v{a)  m{a)  da 


(24) 


S’nce  the  size,  the  fall  velocity  and  the  mass  of 
snow  particles  vary  from  one  snowfall  to  another, 
a  range  of  values  for  k  is  expected,  with  lower  k 
values  being  associated  with  heavy,  fast-falling 
particles. 

Numerous  empirical  studies  I  o  determine  val¬ 
ues  for  k  have  b°en  conducted  by  simultaneously 
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(26) 


measuring  the  extinction  coefficient  Bc«  and  Rs. 
Mason  (1978)  and  Seagraves  (1986)  reviewed  these 
empirical  studies  and  found  the  results  to  .vary 
widely.  A  summary  of  previously  observed  B„x,to 
Rs  relationships  is  presented  in  Table  3.  The  dis¬ 
crepancies  in  these  results  are  probably  the  result 
of  different  snowfall  conditions,  the  different  ge¬ 
ometries  of  the  transmission  systems  and  discrep¬ 
ancies  in  the  snow  precipitation  rate  measure¬ 
ments. 

Table  3.  Relationships  between  measured 
extinction  coefficient  and  snow  precipita- 
ion  rate. 

Bext  =3.2  •  R$  91  Poljakova  and  Tretjakov  (1 960) 

Bext  =4.1  Rs  Lillesaeter(1965) 

Bext=2.4'  Rs  Warner  and  Gunn  (1969) 

Bexl  =3.1  ■  Rs  O’Brien  (1970) 

Bext=1.9  Rs  Nishitsuji  and  Matsumoto  (1971) 

Bex{=2.5-  Rg77  Muer  ch  and  Brown  (1977) 

The  precipitation  rates  for  most  of  the  above- 
mentioned  studies  were  measured  using  mechani¬ 
cal  snow  gauges.  The  poor  time  resolut:on  of  these 
gauges  required  the  data  to  be  averaged  overlong 
periods.  During  the  SNOW  tests,  the  value  fork 
was  reinvestigated  using  electronic-balance  snow 
gauges.  Linear  regression  results  yielded  several 
£c>l  to  R,  relationships,  which  generally  ranged 
from 


Bext  —  1*0  ■ 

Rs 

(25a) 

Bext  =  7.C  • 

Rs 

(25b) 

These  results  are  consistent  with  the  previous 
empirical  findings  (Table  3).  The  slightly  wider 
range  of  k  values  observed  during  the  SNOW  tests 
may  have  been  caused  by  the  shorter  averaging 
time  (1  minute)  used  to  determine  the  precipita¬ 
tion  rate.  The  importance  of  snow  particle  type  on 
the  Bex, to  Rs  relationship  is  evid^pt  by  the  range  of 
values  reported  for  k.  Since  the  Bext  to  Rs  relation¬ 
ship  varies  widely,  the  precipitation  rate  alone  is 
not  a  suitable  parameter  for  use  in  the  evaluation 
of  electro-optical  system  performance. 

Mass  concentration  and  extinction 
The  snow  mass  concentration  is  the  snow  mass 
per  unit  volume  of  air.  It  can  be  approximated  by 


Ms  =  n(a )  in  (a )  da  . 

Jo 

The  relationship  between  the  extinction  coeffi¬ 
cient  and  the  snow  mass  concentration  is 

2na2n(a)da 

Bext  =  - Ms.  (27) 

n(a )  m(a )  da 

Jo 

This  expression  is  similar  to  eq  24;  however,  the 
uncertainty  due  to  variations  in  v(a)  is  eliminated. 
Therefore,  the  extinction  coefficient  can  be  more 
accurately  estimated  from  the  measurement  of 
mass  concentration  than  it  can  be  from  the  snow 
precipitation  rate  measurement.  This  is  illustrated 
in  Figure  20  which  shows  the  relationships  be¬ 
tween  the  extinction  coefficient  and  the  two  snow¬ 
fall  parameters.  The  plot  of  theBex,tc  Ms  relation¬ 
ship  clearly  shows  less  scatter  than  the  Be*:  to  Rk 
plot. 

A 

Empirical  Bcxl  to  Ms relationships  have  been  well 
documented  during  the  SNOW  tests.  Lacombe 
and  O'Brien  (1982)  used  a  linear  regression  and 
found  good  correlation  between  the  visible  and 
infrared  extinction  coefficients  and  the  snow  mass 
concentration.  Their  results  are  shown  in  Table  4. 
Lacombe  et  al.  (1983)  segregated  visible  extinction 
and  snow  mass  concentration  data  according  to 
the  snowjype  and  observed  improved  correlation 
between  B^and  M$.  For  all  types  of  snow,  a  power 
fit  curve  to  the  data  yielded  the  following  relation¬ 
ship 

Bext  -  9.7-  Mg  '88.  (28) 

Hutt  et al.  (1986)  conducted  an  independent  study 
of  the  Be*,  to  Ms  relationship.  A  linear  regression 
analysis  of  their  data  yielded  the  following  rela¬ 
tionship: 

Bext  =k  -  Ms  <29> 

where  k  ranged  from  12.5  to  16.7. 

Snow  mass  extinction  coefficient 
A  parameter  that  is  commonly  used  to  charac¬ 
terize  an  atmospheric  obscurant  is  the  mass  extinc¬ 
tion  coefficient  a  ,  which  is  a  measure  of  how  ef- 

m 
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Figure  20.  Comparison  of  the  Bcj(  to  Ms  and  Ba(  to  Rs  relationships  observed  during  a  snowfall. 


Table  4.  Relationship  between 
measured  extinction  coefficient 
and  .mow  mass  concentration  (af¬ 
ter  Lacombe  and  O'Brien  1982). 


Wavelength 

Bgxf  =k' 

(tun) 

(knr') 

0.55 

Bext  =13.41  Mg 

1.06 

Bext  =14. 16  /vtg 

3.0 

Bext  =19-21 -Ms 

10.37 

Bext  =19.05-  Ms 

Figure  21.  Distribution  of  the  snow  mass 
extinction  coefficient  for  visible  wavelength 
obtained  from  simultaneous  measurements 
of  snow  mass  concentration  and  visible 
transmittance. 


fectively  a  fixed  mass  of  a  particle  scatters  and  ab¬ 
sorbs  the  radiation  incident  upon  it.  The  relation¬ 
ship  between  the  theoretical  transmittance  and  the 
mass  extinction  coefficient  is  expressed  as 


used  eq32  to  determine  the  frequency  distribution 
of  am  that  occurred  during  the  SNOW  tests.  The 
dis  tribution  shown  in  Figure  21  indicates  that  0.03 
m2/g  is  an  appropriate  value  for  am. 


T  =  exp(-am  •  Ms  -  L) 


(30) 


Equation  30  can  be  rearranged  to 


am  =■- 


-In  T 


Ms-L 


(31) 


to  show  that  a  can  be  determined  from  simulta- 

m 

neous  measurements  of  transmittance  and  mass 
concentration.  To  account  for  the  forward-scat¬ 
tered  light  that  enters  the  detector  field  of  view 
during  transmission  measurements,  eq  31  has  to 
be  modified  in  accordance  with  eq  14  and  19  to 


ctm  = 


-InT 


MSL(1-0.5F) 


(32) 


A 

where  Tis  the  measured  transmittance.  Koh  (1986) 


Multiple  scatteringeffect 
As  the  optical  depth  (Btxt  L)  of  a  snow-filled  at¬ 
mosphere  increases,  the  effect  of  multiple  scatter¬ 
ing  eventually  becomes  evident  (nonexponential 
attenuation  of  the  transmitted  energy  is  observed). 
A  two-stream  approximation  to  radiative  transfer 
can  be  used  to  derive  a  simple  analytical  expres¬ 
sion  that  describes  the  effects  of  multiple  scatter¬ 
ing  on  the  transmitted  radiation  as  a  function  of 
the  optical  depth,  the  single  scattering  albedo,  and 
moments  of  the  phase  function.  Bohren  (1987) 
derives  a  simple  analytical  expression  for  the  trans¬ 
mitted  radiation  by  assuming  that  the  light  scat¬ 
tered  by  the  particles  travels  in  two  directions 
only,  exactly  forward  and  exactly  backward.  The 
transmitted  radiation  in  the  two-stream  approxi¬ 
mation  is  expressed  as 
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where  the  asymmetry  parameter  g  is  the  mean 
cosine  of  the  scattering  angle. 

The  transmitted  radiation  is  composed  of 
scattered  Tsand  unscattered  Tucomponents,  which 
can  be  expressed  as 


Ts=  T-Tu. 


(34) 


where 


Tu  =  exp(-Bext  •  L)  ^  (35) 

The  measured  transmittance  T  is  the  sum  of  the 
scattered  and  unscattered  radiation  that  enters  the 
field  of  view  of  the  detector.  The  measured  trans¬ 
mittance  in  a  multiple  scattering  medium  can  be 
expressed  as 

T=TU  +  f  ■  Ts  (36) 

where  /  is  equal  to  the  fraction  of  the  scattered 
radiation  in  the  forward  direction  that  is  measured 
by  a  detector.  The  value  of /is  dependent  on  the 
beam  geometry,  the  optical  depth  and  the  snow 
asymmetry  parameter.  Figure  22  is  a  comparison 
of  visible  transmission  data  obtained  during  the 
SNOW  tests  with  those  calculated  using  eq  36.  A 
more  detailed  derivation  of  eq  36  is  given  by  Koh 
(1989). 


Visibility  and  extinction 
Snow  alters  the  contrast  of  the  atmosphere  so 
that  the  visual  range  Vr decreases  as  the  snow  load 
in  the  atmosphere  increases.  Visual  range  is  de¬ 
fined  as  (Middleton  1952) 

Vr  =  JL_in(C/e)  (37) 

Bext 

where  C  is  the  inherent  contrast  of  the  target 
against  the  background  and  e  is  the  observer's 
contrast  threshold.  A  widely  used  formula  for 
characterizing  the  effect  of  the  atmosphere  on  the 
performance  of  electro-optical  devices  is  the  Ko¬ 
schmieder's  relationship  for  the  meteorological 
range  Vm.  This  is  obtained  by  solving  eq  37  for  a 
blade  target  against  a  horizon  sky  (C  =  1)  and 
assuming  that  e  =  0.02.  Equation  37  then  becomes 


Vr  =  Vm  =  Mi2 

Bext 


where  Vm  is  defined  as  the  meteorological  range. 
The  meteorological  range  is  often  assumed  to  be 


Figure 22.  Comparison  of  thevisible  transmis¬ 
sion  data  with  those  calculated  using  the  two- 
stream  approximation  to  radiative  transfer. 


the  visibility  (how  far  a  person  can  see),  so  that  a 
simple  inverse  relationship  between  extinction 
coefficient  and  visibility  is  assumed.  This  is  of 
practical  interest  since  the  extinction  coefficient 
may  be  estimated  from  how  far  an  observer  can  see 
during  a  snowstorm. 

During  the  SNOW  tests,  experiments  were  con¬ 
ducted  with  trained  observers  to  measure  the  visi¬ 
bility  of  predetermined  targets  during  snowfall 
(Lacombe  and  Petzko  1987).  These  visibility  meas¬ 
urements  were  converted  to  extinction  coefficients 
using  the  Koschmieder's  relationship,  and  then 
compared  with  the  extinction  coefficients  meas¬ 
ured  with  a  transmissometer.  The  results  shown  in 
Figure  23  suggest  that  large  differences  can  occur. 

Although  the  simple  inverse  relationship  as¬ 
sumed  between  the  visiblity  and  the  extinction 
coeffident  makes  the  Koschmieder's  relationship 
an  attractive  tool,  its  practical  application  may  be 
limited.  The  value  of  0.02  for  the  contrast  threshold 
is  controversial  and  large  discrepancies  appear  in 
the  literature.  The  contrast  threshold  is  dependent 
on  the  angle  subtended  by  the  target,  the  sharp¬ 
ness  of  the  boundary  between  the  target  and  the 
background,  and  the  level  of  illumination.  Even  if 
these  factors  are  held  constant,  there  are  differ¬ 
ences  among  observers.  In  addition,  the  ideal  black 
target  against  the  horizon  may  not  be  present 
when  the  visibility  measurements  are  made.  Fi¬ 
nally,  the  use  of  the  Koschmieder's  relationship  is 
limited  to  daytime.  At  night,  a  light  source  of 
known  intensity  is  required  to  make  visibility 
measurements. 
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more  suited  for  military  applications  since  such 
instruments  are  portable  and  easy  to  deploy. 

Measurementof  extinction  by  light-scattering  is 
based  on  the  principle  that  the  intensity  of  the  scat¬ 
tered  light  at  certain  angles  is  linearly  related  to  the 
extinction  coefficient.  Mathematically,  this  rela¬ 
tionship  is  expressed  as 


>ext 


=  constant.  (39) 


,2* 


62 


2  >9(0)  sin9  d0  d<|> 


Figure  23..  Comparison  of  the  extinction 
coefficients  estimated  from  visibility  meas¬ 
urements  with  those  measured  with  a 
transmissometer  ( after  Lacombeand  Petko 
1987). 

Direct  measurement  of  extinction 
Transmittance  measurement  is  the  preferred 
technique  for  determining  the  extinction  coeffi¬ 
cient  through  falling  snow..  However,  such  meas¬ 
urement  requires  a  path  length  of  several  hundred 
meters  to  provide  accurate  results  over  a  broad 
range  of  snowfall  conditions.  Hiis  makes  it  im¬ 
practical  for  military  application,  where  a  simple 
technique  to  directly  measure  the  snow  extinction 
coefficient  in  an  area  where  an  electro-optical  sys¬ 
tem  is  to  be  deployed  is  of  practical  interest.  Ex¬ 
tinction  coefficient  measurement  devices  that  use 
a  short  path  light-scattering  technique  may  be 


Since  detectors  for  measuring  light  possess  finite 
acceptance  angles,  the  angular  scattering  coeffi¬ 
cient  p(0)  is  measured  over  some  range  0,  to  02. 
Theoretical  and  experimental  studies  by  Win- 
stanely  and  Adams  (1975)  have  shown  that  a  for¬ 
ward-scatter  meter  can  be  used  to  measure  extinc¬ 
tion  coefficient  in  fog  using  the  relationship  ex¬ 
pressed  in  eq  39. 

Experiments  by  Muench  and  Brown  (1977)  and 
Koh  (1987)  have  shown  that  forward-scatter  me¬ 
ters  designed  to  operate  in  fog  can  also  be  used  for . 
measuring  extinction  coefficient  in  felling  snow.  A 
forward-scatter  meter  shown  in  Figure  24,  which 
measures  the  light  scattered  between  27  and  47*, 
was  used  to  measure  extinction  in  falling  snow, 
and  its  results  were  compared  with  the  results 
obtained  using  transmission  technque.  Figure  25 
indicates  that  a  foward-scatter  meter,  once  cali¬ 
brated  for  snow,  can  be  used  to  predict  visible  and 
infrared  transmission  through  falling  snow. 


Figure  24.  HSS  forward-scatter  meter  used  for  measuring  extinction  in  falling 
snow. 
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Figure  25.  Comparison  of  the  measured  transmittance  in  the  visible  and  infrared 
wavelengths  with  those  predicted  using  a  fonoard-scatter  meter  designed  for  oper¬ 
ation  in  fog.  A  different  calibration  is  required  for  snow.. 


CONCLUSIONS 

The  effect  of  falling  snow  on  the  performance  of 
an  electro-optical  system  depends  on  the  particu¬ 
lar  system  and  its  intended  function.  Snow  in  the 
atmosphere  scatters  and  absorbs  optical  radiation 
and  therefore  degrades  the  performance  of  many 
sensors  operating  in  a  winter  environment.  The 
SNOW  tests  were  conducted  to  investigate  the 
effect  of  snowfall  on  the  optical  properties  of  the 
propagation  path.  Numerous  transmittance  meas¬ 
urements  in  the  visible  and  infrared  regions  were 
made  in  conjunction  with  detailed  measurements 
to  characterize  the  falling  snow.  As  a  result,  new 
empirical  and  theoretical  models  to  predict  the 
attenuation  of  optical  radiation  in  falling  snow 
and  improved  techniques  to  characterize  the  prop¬ 
erties  of  falling  snow  have  been  developed. 

More  studies  are  required  to  develop  a  satisfac¬ 
tory  analytical  expression  to  describe  the  snow 
size  distribution.  Our  inability  to  model  the  com¬ 
plex  mechanism  of  snowfall  formation  and  to  accu¬ 
rately  measure  the  size  of  snow  particles  has  lim¬ 
ited  snow  size  distribution  studies.  Information 
on  the  phase  function  of  snow  particles  is  also 
limited.  Experiments  are  required  to  measure  the 
angular  light-scattering  properties  of  snow  to  more 
effectively  address  the  effects  of  single  and  mul¬ 
tiple  scattering  on  the  energy  measured  by  any  op¬ 
tical  system. 

Although  much  effort  has  been  made  to  study 
the  physical  and  optical  properties  of  fallingsnow, 
it  is  necessary  in  the  future  to  measure  the  effects 
of  the  snow-filled  path  in  terms  of  the  parameters 
that  are  useful  in  the  analysis  of  electro-optical 
systems.  The  extinction  coefficient  is  a  general 


measurement  of  attenuation  that  can  be  applied 
uniformly  to  all  electro-optical  systems.  However, 
system  performance  predictions  may  require  para¬ 
meters  beside  the  extinction  coefficient.  Future 
study  of  the  optical  properties  of  the  snow-filled 
atmosphere  that  address  such  parameters  as  the 
modulation  transfer  function,  the  point  spread 
function,  and  time  dispersion  of  the  snow-filled 
atmosphere  will  be  useful. 
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